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Abstract 

We examine the problem of the dynamics of interfaces in a one-dimensional 
space-time discrete dynamical system. Two different regimes are studied: the non- 
propagating and the propagating one. In the first case, after proving the existence 
of such solutions, we show how they can be described using Taylor expansions. 
The second situation deals with the assumption of a travelling wave to follow the 
kink propagation. Then a comparison with the corresponding continuous model is 
proposed. We find that these methods are useful in simple dynamical situations but 
their application to complex dynamical behaviour is not yet understood. 
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1 Introduction 



Coupled Map Lattices (CML) have been proposed as simplest models of extended 
dynamical systems and they serve now as a paradigm in the analysis of complex 
nonlinear phenomena. CML is a discrete space, discrete time, continuous state 
dynamical system that is deterministic [Q || ||. These models are easy to handle 
numerically but the mathematical analysis is not so simple. Also they often display 
the essential features of physical systems (for a review, see for example [§]). The 
numerical studies show that CML give rise to a rich phenomenology including a 
wide variety of both spatial and temporal periodic structures, intermittency, chaos, 
domain walls, kink dynamics, etc. || ||. Some of these phenomena were also 
investigated from a rigourous mathematical point of view [§] || . 

The problem that we are dealing with is the study of particular localized coherent 
structures in CML, the interfaces between two different phases, the so-called kinks 
or fronts. Depending on the nature of the phases, these interfaces can emerge in a 
simple or rather complex dynamical context. In the former case, kinks appear when 



both (stationary) phases are in competition as is the case in alloy solidification (1C 



or in the propagation of chemical waves [IjJ . Interfaces in such systems can be either 



propagating or non-propagating fronts depending on the stability of both phases and 
on the value of the control parameter, here taken as the diffusion coefficient. The 
dynamical behaviour is said to be simple because of the low temporal period and 
the (quasi-) homogeneity of the different phases, namely the domains. On the other 
hand, in spatiotemporal intermittency, phases are more complex solutions (in the 
sense of periodicity) and they have been called the natural wavelengths. Kinks in 
this case play again a central role as their motions destabilize these wavelengths and 
thus lead to the well-known bursts of turbulence inherent in intermittency || fll2| . 

The goal of this article is to give a rigourous description of the dynamics of the 
kinks that may occur in CML models. The advantage of CML is based on the use of 
local maps with well-controled dynamical properties. Tuning the parameters in such 
systems enables the observation and study of various regimes of the front dynamics. 
Moreover it gives a discrete alternative point of view to the more conventional tool 
of partial differential equations (PDE) in which the dynamics of localized structures, 
at least for some PDE, begins to be well understood [|13| flU}] . 

We shall compare the behaviours of these interfaces in the discrete space-time 
dynamical system and the corresponding continuous one. After introducing the 
model and patterns under consideration, we concentrate on the situation of non- 
propagating solutions in Section 3. We show the existence and unicity (up to trans- 
lations in the lattice) of such kinks in the small coupling range. Its expression is 
then evaluated using Taylor expansions. In a second step, the propagating fronts are 
analyzed with the assumption of being a travelling wave as was already done in PDE 
|T5| . In this case there are important differences between discretized and continuous 
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models concerning the dynamics of interfaces and we address this issue in Section 
4. Our assumptions are tested numerically. Finally an extension of the preceding 
results to upper temporal period is presented and we end with some concluding 
remarks. 



2 Definitions 

In this section we begin the description the model using the general definition given 



in 
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Let ^ be the (infinite) unidimensional lattice. At each site i G 'E the local phase 
space is Xi (here Xi is chosen as either the interval [0, 1] or [—1, 1]). The phase space 
of our dynamical system is the direct product M = Xi. Consequently a point 

x G M can be written x = (xi) ie ^. 

CML is a mapping $ e of M into itself depending on the coupling parameter 
e. As we are dealing with a model of reaction-diffusion systems, $ e splits into the 
composition of two mappings 

$ e = G e oF (2.1) 

where (Fx)i = f{xi) is the reaction, that is to say the action of a local nonlinear 
mapping / : Xi — > Xi and (G e x)i = g e ,i( x ) * s the interaction. In our case this 
interaction is chosen to be representative of the diffusion in the simplest way which 
reads 

g e ,i(x) = Xi + ~(x i+1 + - 2xi) 

where the coupling parameter e G [0,1]. Finally the new state of the CML at time 
t + 1 is given by the following convex linear combination (of the state at time t): 

xf l = (<S> € x% = (1 - e)/(^) + \ (f(xl +1 ) + f{x\_ x )) (2.2) 

which implies the location of x t+1 in M as soon as x l G M. 

Let {X 1 , X 2 , • ■ ■ , X p } denote either p different fixed points or a period-p orbit of /. 



Definition: A left domain (respectively a right domain) is a sequence x = {a^}^-^ 

(resp. x = {xi} % ~^) such that there is an X p such that 

\x\ - X p \ < a for all % < Nt (resp. \x\ - X p \ < a for all % > N 2 ) 

where a is a small positive parameter to be fixed later. 

The different values X p are referred to as phases of the domain. A kink is the 
interface between a left and a right domain with different phases (fig. 1). A rigorous 
definition of these patterns can be given with a = 0, but this implies the use of 
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local maps that are constant over an open subset of X{ and this constant is a point of 
the considered periodic orbit. In our case, we deal with more general maps that do 
not have these peculiar features. Thus, it is not possible to find kink-like solutions 
with domain values (strictly) equal to one of the X n . In other words domains do not 
(really) exist when the lattice reveals the presence of a kink and when e is non-zero. 
This explains the above definition where a can be thought of as a control parameter 
for separating the interface from the coherent domains. The smaller is a, the longer 
is the interface. 

3 Non-propagating kinks 

We address in this section the problem of non-propagating interfaces. After proving 
the existence of such kinks in the small coupling range, we show how the use of 
Taylor expansions allows us to describe these solutions and induces the reduction of 
the dynamics to a finite dimensional space. 

For the sake of simplicity, the map / is assumed to be fc-times continuously 
differentiable with bounded derivatives and also to have two stable fixed points X\ 
and X 2 . Thus the domains under consideration are fixed structures for the CML 
dynamics and the kinks are inspected as fixed points. 

We are looking for the solutions of <& e (x) = x with the boundary conditions 
lim Xi = X\ and limxj = X 2 . When e = there exists a trivial solution x* and it 

i^r— oo i— >oo 

is easy to see that x* is invariant under translations on the lattice. This solution is 
written 

x* = (•■■, X\X\X\X\---) 
Now we endow the phase space M with the sup norm: ||x|| = sup|xj|. One can check 

that (M, ||.||) is a closed ball of a Banach spaceQ. Moreover the Frechet derivative 
of $ at x*, £>$ ( 

x ) , is ci bijective operator from M into R since it is diagonal 
and none of its elements vanishes. Finally the operators <3> e and -D$o axe bounded 
and continuous at x*. 

Thus by the Implicit Function Theorem, we obtain the existence and unicity of the 
kink-like fixed point in the neighborhood of e = 0, that is to say in a certain interval 
of strength of space interactions [0, e'[. We may note that, in the more general case 
of temporally periodic domains of period p, this theorem can also be invoked using 
($ e ) p instead of $ e to prove the existence and unicity of the kink-like periodic-p 
orbit. 

Since we know the existence of the interface when e is small, we can now describe 
it. A first step in this problem is to follow the kink genesis from a perturbative point 

*We will keep the notation M for this Banach space. 
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of view. The system is initialized in the following state: 



Vz < s 
Vi > s 



x 



X 1 
X 2 



(3.1) 



where s is fixed in Z . (3.1) is the so-called step. 

This pattern is submitted to the dynamics and the successive states are computed 
using expansions (close to X 1 if i < s, to X 2 otherwise) in terms of e powers. When 
these expansions are limited to the k th power, it can be proved that x 1 (t > k) is 
given by the system: 



X j 

x\ 

x i 
1 1 

X f ; 



X 1 
X 1 



o e" 



X 1 ^ a itl (t)6 l + o(e k 

l=s-i+l 

X 2 + £ a iAty + o{e k ) 

l=i — s 

X 2 + o(e k ) 
X 2 



-t + s 
-k + s 



k + s 
t + s 



i < 

< i < 

< i < 

< i < 

< i < 

< i 



-t + s 
-k + s 



k + s 
t + s 



(3.2) 



where o(e k ) stands for the remainders of the expansion. In this expression, a,ij(t) 
is the I th power coefficient for the polynomial determining x\. This general sys- 
tem manifests itself for any local map and the diffusive features only enter in the 
coefficients. It reveals a kink where both domains are defined by a = o(e k ). 

Such a solution can be understood in the following way: in the small coupling range 
the perturbative effects (from the local map point of view) induced by the presence 
of the interface are rapidly damped in space. This explains the limited number of 
polynomials appearing in the characterization of the structure. 

From the system (3.2) we deduce the interface length expressed in terms of the 
number of intersites: Ld(k) = 2k. In order to display this kink length as a function 
of the coupling parameter, one has to use the relation a ~ e k . We obtain obviously 



2 Log a 
Loge 



We describe here a first comparison with continuous systems. From a PDE model 



the interface length reads [14| 



Thus both lengths vanish when e goes to zero. However their comparison tells us 
that L c decreases more rapidly than L^. There then appears a difference between 
the models, starting without diffusion and increasing slowly the coupling parameter, 
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according to both definitions of an interface, we get a longer kink in a continuous 
model than in a discrete one. 

Furthermore the structure of (3.2) is invariant under the CML dynamics, that is to 
say only the e polynomial coefficients evolve in time. In fact, for any given k and s, 
substituing the corresponding (3.2) expression in (2.2) using the Taylor expansion 
(until the k th power) for / proves the invariance. 

The kink pattern is therefore given by a set of polynomials with time-dependent 
coefficients. So when studying such structures, the CML dynamics can thus be re- 
placed by a map defined on these coefficients. This substitution allows the reduction 
from an infinite number of relations (in the latter dynamics) to just a few (in the 
coefficient map). In other words the stability problem is now located in a finite di- 
mensional space which is an issue to the closure problem in such unbounded media. 



The coefficient map analysis is made here for s = and k = 2 since it is the simplest 
case and since the other cases can be treated in the same way. When the expansions 
are limited to the second power, there are six non-vanishing coefficients for which 
the dynamic F is (defined in R 6 ) 



o_ 1>2 (t + l) = (a_ 1>2 (t) + ^/'pT 1 ) 

ao,i(* + l) = X + ao^f'iX 1 ) 

a 0>2 (t + l) = ( ao>a (t)-a 0> i(t))/ / (X 1 )+a 1>1 (t)£^ + a g >1 (t) £ ^ 

a 1>2 (t +1) = (a 1;2 (i) - a 1;1 (t))f(X 2 ) + a 0jl (t)^p + aj^^p 

a 1A (t+l) = -X + a 1A (t)f'(X 2 ) 

a 2 , 2 (t + l) = (a 2 , 2 (t) + a -^)f'(X 2 ) 



(3.3) 



X 2 — X 1 

where X = and /', /" denote successives derivatives with respect to x. 

Expression (3.3) is deduced from the introduction of the system (3.2) in the CML 
dynamics. The iterations for a_i )2 and ao,i (respectively a 2)2 and a^i) are in this 
particular case the relations linking the coefficients. 

Before looking at the map F, we make a comment on the influence of / orbits 
on the kink dynamics. If, instead of choosing a local map with two fixed points, we 
work with a period-2 orbit {X 1 , X 2 }, then the interface would flip at every iteration 
as the domains' values alternate. The kink structure would therefore be given by 
the system (3.2) for even times and by a symmetrical (with respect to the index s) 
expression for odd times. Such a solution describes a period-2 orbit for the lattice 
dynamics. But as the coefficients' expressions in this case are quite complicated we 
refrain from presenting them here and we will just mention the final results. 

Notice that these expressions are s-independent because of the translational invari- 
ance in the lattice frame. 



5 



The analysis of the map F shows the existence of a fixed point called A°° for 
which the coordinates are 



a°° 



2(1 -/'(X*)) 2 
X 



*o,i 



*0,2 



1 - f (Xi) 

-X ( /'(X 1 ) f(X 2 ) xfix 1 ] 



+ 



1- /'(Xi) VI -/'(X 1 ) 2(1 -/'(X 2 )) 2(1 -f(X*)) 2 

* (JV^l_ , fix 1 ) x/"(x 2 ) 



l-/'(X 2 )Vl-/'(X 2 ) 2(1 -/'(Xi)) 2(1 -/'(X 2 )) 2 
X 



"Li ~ i _ /'(x 2 ) 



"2,2 



2(1 -r(x*)y 



One can find „4.°° when looking for a solution of A = F(A), A e R 6 . 
We now compute the jacobian matrix associated to A°° as its spectrum determines 
the stability of .4°°. When solving the v4°°-spectral problem we found three eigen- 
values equal to /'(X 1 ) andthree others equal to /'(X 2 ). Since X 1 and X 2 are two 
stable fixed points, ^4°° is stable. Moreover the map F is strictly contracting so this 
fixed point is unique. Thus the expression (3.2) defines the unique attractor of the 
CML dynamics for the considered expansions. 

To conclude, the kink stability is only due to local map stability; once the domains 
are installed, there is a unique stable pattern possible for the kink as much as small 
diffusive interactions are taken into account. This is important because it proves 
that no special preparation is needed to create and stabilize such interfaces; they 
arise in a natural way, no matter wath the initial condition. One can see that this 
technique not only confirms the results given by the Implicit Function Theorem 
used at the begining of this section, but also gives an explicit expression for the 
patterns under consideration. The validity domain of these solutions will be made 
more precise later on with a localization argument. 

A similar calculus, pursued in the case k — 3, generates a coefficient dynamical 
system in R 12 and leads to the same conclusion. This technique can also be used 
to examine the case of greater periodicity as said before. One has to compute 
expansions in different phases of the dynamics. Other coefficient maps are found. 
We summarize here some of the results. When / has a periodic orbit of period two 
(for k = 2 and k — 3) the jacobian matrix associated to the CML periodic orbit is 
shown to have only one eigenvalue, namely /'(X 1 )./'(X 2 ). In the case k = 2 this 
eigenvalue has multiplicity 6 and 12 when k — 3. The map is also contracting which 
allows in these cases the same conclusions as above. An extension of this result to 
any finite temporal period needs more calculus and is still in progress. 
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4 Propagating fronts 



We now consider a particular aspect of the kink dynamics, namely its displacement 
along the lattice. As already stated, this phenomenom appears as one of the keys 
of the spatio-temporal intermittency in CML. As in the fixed interface situation 
we begin by investigating the case where / has two stable fixed points, then we 
generalize to greater temporal periods. 

A first step in this direction is made using computer simulations. They suggest 
that these interfaces are moving without being deformed and with a constant velocity 
when the local behaviour is quite simple. Starting from this statement, we assume 
the existence of such a kink-like travelling wave in CML. If that is the case, the 
lattice state at the date t is given by 

x\ = h{% -vt) Vte NVie % (4.1) 

where we suppose h G C 1 (l?, /) and v is the travelling velocity. For the present 
purpose, assumptions for the local map are kept identical with those made in Section 
3. The continuity implies the existence of a third fixed point X u (X 1 < X u < X 2 ) 
which is unstable [[H]]. Likewise we impose the boundary conditions for h in the 



co-moving frame z = i — vt: 

lim h(z) = X 1 , lim h(z) = X 2 

z— »— oo 2^+00 

We call z u the point such that: X u = h(z u ). 

If a slow motion (as seen in simulations) is considered, that is if v <C 1, in one 
time step the site values Xi stay close to each other. Consequently we write in the 
first order approximation 



h(z — v ) = h(z) — v^— 
dz 



(4.2) 

2=1 — Vt 



In a second step the two right-hand sides of the relations (2.2) and (4.2) are com- 
pared. The expression obtained by generalization to all real values of the variable z 
is 

h{z) - vh'(z) = f[h(z)} + ^Af[h(z)} (4.3) 
Eqn. (4.3) is now intergrated over R. From the h hypothesis we deduce 

/+oo 
h'{z) dz = X 2 - X 1 
-oo 
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and: 

/+oo 
Af[h{z)\ dz = 
-oo 



Finally the moving speed is given by: 

1 



X 2 -X 1 



(h(z) - f[h(z)}) dz (4.4) 



This relation indicates that the kink displacement is a consequence of a lack of sym- 
metry in the interface shape and its image under /. We also notice from (4.4) that 
a necessary condition for v to vanish is that the functions / and h are antisym- 
metric functions. A similar calculus performed in a continuous space-time model of 
chemical reactions had lead to similar conclusions IT5 . 



Following these comments, this inhomogeneous solution is now studied for a 
peculiar local map, say, the cubic map f(x) — v x (1 — x 2 ) (1 < v < 2). As this 
reaction is antisymmetric, the kink is expected also to be antisymmetric and thus 
to be a stationary solution. In order to compute an analytic expression for this 
pattern, we use a continuous description of our model assuming kinks in CML can 
be described by the following equation (See e.g. [|18|1 ) 



e d 2 h(z) 



h(z) = f(h(z)) + -^y- (4.5) 

One can find the solutions of (4.5) using an analogy with the equation of motion of 
a particle of mass | at position h, in time z evolving in the potential 

V(h) = jf (/(*») - h)dh = ^h 2 - v -tf 

Investigating a first integral of motion (the constant is given by the above boundary 
conditions) leads to the kink profile in the stationary case 




h(z) = \ tanh W z) (4.6) 



A plot of this function verifies the validity of the approximation (4.5) (fig. 2). One 
can obtain a better fit with the discrete kink when considering instead of (4.5) the 
equation 

r / \ f / 1 / \ \ e d 2 f(h(z)) 

The mechanical analogy is also applicable in this model but no simple analytic 
expression for h was found. 
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In any case we take (4.6) as a starting point for understanding the non-symmetric 
case and the single map is chosen to be the perturbed cubic map f(x) — ux(l — 
x 2 ) + p (p <C 1). We suppose that the effect of the perturbation is to induce a 
displacement of the interface. So we make again the ansatz of the travelling wave 
(which has been proved for a certain class of partial differential equations pJ|): 




h(i t t) = Cp+\ tanh W (i-p.v.t)) (C e R) (4.7) 



Eqn. (4.7) is substituted in (4.5) also supposing a small speed. For the first order 
term this yields 

3 Jev 

v = 4.8 
2 v — 1 

It is also possible to compute the constant C but, since it plays a secondary role in 
this solution, we do not give its expression here. More interesting is the comparison 
among the expressions (4.4) and (4.8) and the measured speed from numerical simu- 
lations (fig. 3). The gaps between the experimental data for the velocity are due to 
the finite size of the lattice used in the numerical experiment. Notice that there is an 
acceptable fit for large e, while we can see a difference between the continuous model 
and the CML behaviour when the diffusion plays a small role. This is because the 
approximation made in (4.5) is no longer valid when e decreases, as discrete effects 
become more and more important. In other words, the discrete laplacian opera- 
tor cannot be considered as a diffusion operator in the classical sense in the small 
coupling range as its variations along the front are much more important and, in a 
certain sense, less regular than those of the corresponding continuous operator. Fur- 
thermore this figure also shows that the assumption (4.1) concerning the existence 
of a travelling wave in CML is in good agreement with the experimental results for 
any strength of space interactions. This is confirmed by a plot of the velocity distri- 
bution when the kink propagates in space and by the expansion method in the fixed 
case. Indeed, it is proved that computing the velocity v using (4.4) and the interface 
shape given by system (3.2) with asymptotic coefficients always gives a vanishing 
value. However the proof of the existence of such a travelling wave solution in CML 
is still an open problem. 

We may understand in a more detailed way the difference between CML and 
a continuous model. One has to display a discrete argument to explain a non- 
moving interface when the strength of space interactions is small. With the above 
assumptions on the single map, there exists an attracting neighbourhood near each 
fixed point such that / restricted to these open subsets is stricly contracting. We 
shall see that the interface is stationary if each x\ remains located in one of the 
intervals. The following points are defined: Xi n f (resp. Xs up ) such that f'(Xi n f) = 
1 and X 1 < Xi nf < X u (resp. f{X Sup ) = 1 and X u < X Sup < X 2 ). Now we 
consider 0\ —]X x ,Xi n f[ and 2 = }Xs up , X 2 [ as these open attracting subsets: 
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Proposition: There exists an e c such that for all e < e c we have for all t > t 

\/i<N x\ e Ox 
\/i > N x\e0 2 

if this condition is true at t — to and if / is an increasing function on ]X 1 , X 2 [. N 
depends on the initial condition. 

The proof is direct by induction. It is enough to write the different inequalities at 
time t and to submit them to the dynamics using relation (2.2) and the properties: 
f(X In f) < X In f and f(X Sup ) > X Sup - A necessary condition for the location in 
Oi (respectively in 2 ) gives the value: ei = ^rjpT^y (resp. e 2 = ^^H'^gr )- 
Finally we get e c = min{e!, e 2 }. 

In the case of the perturbed cubic map we compute the value e c ~ 0.13 which is in 
good agreement with the experimental data. From this proposition we also obtain an 
effective range of application of the expansion method. This value e c may also serve 
as an upper bound for the e' given by the Implicit Function Theorem. This bound 
is up to now the only mark in the study of the transition from a non-propagating 
to a propagating interface in our model (with a continuous nonlinear map). 

In a more general dynamical context, that is to say when / has a period-]? orbit, 
the analysis can be done in a similar way. In this situation we assume the existence 
of a periodic sequence of travelling waves {hi(z), ■ ■ ■ , h p (z)} describing the CML 
states, each hi(z) being again continuously differentiable. The generalized boundary 
conditions are 

' h+°° = lim hn(i) = X n+1 

i— *+oo 

h-°° = lim h n (i) = X n 
with the periodic condition: X p+1 = X 1 . 

Following the technique employed in the simplest case and assuming that the speed 
is independent of the phase of the solution then yields 



/+oo 
(h n (z) - f[h n (z)}) dz (4.9) 
-oo 



We have tried to test this expression in cases where the period-2 and -4 orbits of the 
logistic map are stable. In both cases the system reveals a non-moving interface for 
any value of the diffusive coupling. This unexpected observation was confirmed by 
the velocities computed from (4.9). This computed velocity vanishes whatever the 
phase n we consider in (4.9). 

We end the study of propagating kinks with a remark on the most complex be- 
haviour, that is to say, an interface between natural wavelengths. In order to obtain 
such a situation, one has to choose a single map with an attracting set having two 
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disjoint intervals exchanging their dynamics P| |TI|. Then the preceding reasoning 
can be pursued with slightly different boundary conditions that are periodic in time, 
and it gives an expression very similar to (4.9). However a comparison with numer- 
ical experiments is quite difficult in this realistic situation. We have noted from 
simulations that the spatial oscillations induce shape modifications of the moving 
structures, the wavelengths destabilize themselves and irregular behaviour is ob- 
served during a transient before the stabilization of a (different) wavelength. Thus 
all the assumptions required of h are no longer valid. For instance, it is easy to 
see that the shape is deformed under iterations and the differentiability is not re- 
spected. We found however some examples of agreement between theoretical and 
numerical results. They are not given here since they do not constitute an effective 
confirmation of the hypothesis. 

5 Summary and discussion 

In this paper we have considered the problem of a peculiar type of inhomogeneous 
solution in CML. First the use of a technique based on expansions to describe 
stationary interfaces has revealed itself as an efficient tool and we concluded that it 
yields to a stable solution in all case. The second goal was to show how a travelling 
wave could properly describe propagating kinks. In order to obtain an analytic 
expression for its shape, a continuous description has proved its usefulness even 
if it is no longer valid in the small coupling parameter range. At first sight one 
can conclude that there is a contradiction between both results, but a localization 
argument was used to link these apparently conflicting assertions. 

We have demonstrated the efficiency of two different tools in the study of the kink 
dynamics in CML. When compared to the case of the PDE, we have seen that the 
effect of the discretization is to shift the critical value for the transition between 
stationary and travelling kinks, from e = in PDE to e c > in CML. In other 
words there appears a vanishing speed in the small coupling range even if / has no 
peculiar symmetry. Our method was shown to generalize to the case of period-p 
orbits. Finally in spatio-temporal intermittency, the validity of such a description 
is still an open question. 
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Figures Captions 

Figure 1: Snapshot of a CML of 30 sites (horizontal axis) after 10000 iterations 
with e = 0.25. The local map is the logistic map f(x) — u, x (1 — x) with \x = 3.2. 
The dashed lines represent domains values, namely X\ and X2. The lattice has been 
initialized with the so-called step (see text). 

Figure 2: Snapshot of a CML of 25 sites (horizontal axis) after 10000 iterations of 
the cubic map, v — 1.4 and e = 0.6. The dashed line is the corresponding continuous 
solution (4.6) 

Figure 3: Plot of the travelling velocity v v.s. the coupling e. The local map is the 
perturbed cubic map with v — 1.3 and p = 0.02. Squares stand for the measured 
velocity obtained with a lattice of 100 sites, the solid line is the one computed (nu- 
merically) by the integration formula (4.4) while the dashed line represents theorical 
velocity obtained in the continuous model, Eq (4.8). The negatives values stand for 
a front propagating toward the left. 
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